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This paper is devoted to some paradoxes occurring in measurements on 
metallic surfaces. They mainly concern the great discrepancies of meas- 
ured values of the electron work function onto monocrystalline surfaces, 
the dependence of this work function on Miller’s indices of walls of crys- 
tal, some misleadings in classic definition of the contact potential and the 
nature of specific adsorption forces which seemingly negate the electro- 
static laws. 

In the paper the calculated values of the work function of noble metal sin- 
gle crystals are presented and discussed with respect to experimental dis- 
crepancies. It is also shown that corrections introduced to work function 
definition by means of Green Function Method (GFM) and Extended Co- 
herent Potential Approximation (ECPA) remove the paradoxes mentioned 
above. The separation of the electrode charge between a bulk and surface 
of the electrode by means of GFM, fully explain the nature of the specific 
adsorption forces. 


1. INTRODUCTION 


Experimental investigations on metallic surfaces have turned out to be sur- 
prising in some cases and they sometimes present big discrepancies in the 
measured values. 

Particularly controversial problem constitutes the significant deviations of 
the measured values of the electron work function of the noble metal (Cu, Ag, 
Au) monocrystallic faces (100), (110) and (111) [1, 2]. Later on these meas- 
urements will be analysed. 

The idea of “an agreement of the calculated values with experiment” com- 
monly used as a test of the relevance of any theoretical model loses in many 
cases its physical meaning due to the lack of an appropriate experimental data 
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in comparison with the numerical calculations. Therefore, any theory which is 
able to produce the numerical results with the highest accuracy constitutes the 
only reasonable source of information of the values interested in. 

In this paper some numerical results of the Green Functions Method (GFM) 
and Extended Coherent Potential Approximation (ECPA) are presented which 
are able to replace the uncertainty of experimental data and, moreover, to ex- 
plain some experimental facts being seemingly in contradiction with the well 
established and recognized laws of physics or chemistry. 


2. PARADOXES IN SOME EXPERIMENTS ON METALLIC SURFACES 


The most controversial experimental results on metallic surfaces concern the 
electron work function and contact potential as well as the so-called the specific 
adsorption of ions on metallic surfaces. 

They can be summarized in the following points: 

1. The great deviations of measured values of the electron work function from 
metal into vacuum (denoted as ġ™ ) are observed for Cu, Ag and Au (100), 
(110) and (111) monocrystals as well as for polycrystalline Cu, presented in 
Table 1 [3-22]. 

Such great discrepancies of the results for the monocrystalline surfaces of 

Cu, Au and Ag seem to be particularly surprising due to the modern tech- 

nology of preparation and cleaning of the samples. 


2. Strong dependence of the measured work function of the single crystal sam- 
ple with different crystallographic faces on Miller’s indices of the faces 
takes place [1, 23]. However, it remains in contradiction with the fact that 
the Fermi level (denoted as Ep) of the measured sample is only one. Tradi- 
tionally, the electron work function is defined as 6 =-E¥ . 


3. The paradox occurs in classical definition of the contact potential applied to 
electrochemical problems. 
The contact potential A@ between two samples of different metals (known 


also in electrochemistry as Volta potential) is defined as the respective dif- 
ference of the work functions of these metals [24, 25] in the form 


Apn = Mel _ Mell — A Mel _ Mell )= (em! I paa) (2.1) 


pe Poze 


where €* stands for the so-called potential of zero charge (pzc) of metal 
Me, i.e. the potential generating the surface charge of metal Q, =0; Q% 
denotes the outer potential of the metallic phase œ in pzc point satisfying 
one of the most fundamental relations of electrochemistry 
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Table 1. Comparison of the work function experimental data for Cu, Ag and Au (100), 
(110), (111) single crystals as well for polycrystalline Cu. »% stands for an average 
value of the data in appropriate column 


4.59% 9” =4.79 4.48" 9” = 4.46 4.90" 9” =4.94 


Polycr. Cu: 4.352"; 4.447"; 4.448", 4.542"; 4.651" o” =4.488 


inkl a ened EE al 
EE E aot bmunmbonaaie ach AoAeteowe bayied ec | 
Ie A Tk el cen mer reese arearen nell 
Dee a Pent = Span tenprarereenn el 
5.00' 9” =5.23 5.13” 9” =5.21 
ME ee ae aoa la o a 


a. Ref. [3] h. Ref. [10] o. Ref. [14] 
b. Ref. [4] i. Ref. [11] p. Ref. [15] 
c. Ref. [5] j. Ref. [17] r. Ref. [12] 
d. Ref. [9] k. Ref. [18] s. Ref. [13] 
e. Ref. [7] l. Ref. [19] t. Ref. [21] 
f. Ref. [6] m. Ref. [20] u. Ref. [22] 
g. Ref. [8] n. Ref. [16] 
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y=" +x" (2.2) 


In eq. (2.2) y* is recognized as an inner potential of the metallic phase a 
(Galvani potential) and y“ stands for a surface potential. 

Formula (2.1) leads really to an apparent paradox. According to the defini- 
tion of classical electrochemistry in pzc point the outer potential of the me- 


tallic phase is equal to zero, i.e. Ø}. =0 [24]. So for two different metals 


Me I and Me II from one side we have Ep” —e" =o" —9"" 40 (be- 


cause of the different work functions of different metals). On the other hand 
leg’ en" J= lop! -pu )=0. 
It is worth to note once more that this paradox is caused by means of the as- 
sumption that in pzc the metallic phase lost its electric charge and then 
Pp =0. 

4. The mechanism of the so-called specific adsorption of ions on metallic 
electrodes. 
This phenomenon was experimentally performed already at the beginning of 
the twenties [26, 27] and from that time it is well known that some anions 
(for instance halides such as CI, Br, I) adsorb themselves on negative 
charged surfaces of the electrodes and some cations (such as Cs*) on posi- 
tively charged ones seemingly against the electrostatic laws. An explanation 
of this unusual phenomena was impossible for many years till the moment 
when the Green Functions technique was applied. 


3. A BRIEF OUTLINE OF THE GREEN FUNCTIONS METHOD AND THE 
EXTENDED COHERENT POTENTIAL APPROXIMATION 


The GFM used in the ECPA leads to numerical calculations of the most im- 
portant quantities of the metallic surfaces such as: occupation numbers of an 
electron, work function into vacuum, surface charge of an electrode, its polari- 
zation by an external voltage and boundary conditions treated as the surface 
potential barrier. Such calculations for noble metals have been carried out in 
earlier papers [1, 2, 28, 29]. For the sake of a better understanding of this paper 
a brief outline of the GFM and ECPA is presented here. 

The Green Functions used in this paper are called the Energetic Green 
Functions (EGF) being the Fourier transformates of the time-dependent re- 
tarded and advanced Green Functions (GF) exhaustively applied in sixties as 
the propagators of the Quantum Field Theory [30-33]. 

Let us now to define the local density of electronic states as the second or- 
der EGF of an electron by means of the formula [28]: 


An application of the Green Functions Method... 237 


G,;(E)=((A,()|8;@)), (3.1) 


where À denotes the quantum state of the electron strictly defined later on and A, 
B are, in general, certain time-dependent operators sometimes used also for t = 0 
as a special case. The operators in eq. (3.1) can be replaced by the chain of the 
appropriate operators giving the higher order GF’s. A concrete form of the GF 
always depends on the Hamiltonian of the problem. 

The Local Density of Electronic States (LDOS) can be expressed as the 
imaginary part of the GF [28] 


p(E.A)=-—ImG,, (E) (3.2) 


So the important achievement of the GFM is that it provides an analytical 
expression for DOS if Im G is assumed to be previously known. The Green 
Function satisfies the equation of motion which in the most general form can be 
written as [30] 


BG,;(E)=——((A, 0,8, ©) +(([A. #71 |; 00), 6.3) 


where the symbols [A, B] and [A, H]- stand for the anticomutator and the co- 
mutator, respectively, and H denotes the Hamiltonian of the considered system. 

In the case of a thin metallic film, consisting of n monoatomic layers, which 
is a quite realistic model for a monocrystalline electrode, the electronic Hamil- 
tonian in the second quantification representation can be used [28] 


À ; - 1 > i + + 
uM AAD 


where A=\,j,,m;o) is the quantum state of an electron in configuration 
(Wannier’s) representation; v =1,....n labels the consecutive number of a layer 
in a film; jy is the position vector of an electron in layer v; m denotes the 
electronic band and o =T (L) is the spin orientation; aj and a, are the crea- 
tion and annihilation operators of an electron, respectively; and finally t,; and 
Izz Stands for hopping and Coulomb or exchange integrals. 

If we apply now the equation of motion (3.3) to the Hamiltonian (3.4), 
making the following identification of operators as A, (0)=at,B,(0)=a,, and 
after decoupling of the fourth order GF appearing in such an equation by means 
of the second order GF [28] one can obtain 
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1 
EG,;(E)= spo +9 K,,G,;(E) (3.5) 
where the matrix elements are taken in the Hartree-Fock approximation 


Ki =t + 2 (Dar -Laiz Mp vp (3.6) 


satisfying the set of secular equations 


F KETY =E (3.7) 


where energy eigenvalues E? and eigenfunctions T}, represent the quantum 


state of the electron in momentum (Bloch’s) space recognized by means of an 
index p. 

In order to simplify the calculation of the GF from eq. (3.5) the ECPA was 
applied [1, 2]. This method consists in the replacement of any real inhomoge- 
neous medium (alloy, melting salt, metallic crystal with spins up and down, 
etc.) by means of an effective homogeneous system, described by an effective 
Hamiltonian, with the so-called coherent potential © pi 

There is a simple relation between the Green Function G,, (E) and coherent 
potential ©, which allows us to substitute ImG,, (E) by Im, [28]. Then the 


final formula of the LDOS can be expressed as 


1 
2 2 1 E-E? \/? 
A AE a |e 3.8) 
p(E,A) gol | l B, | ( 
with the value of a broadening of spectral lines of electrons in the form 
B, = Lil! iMayeyftaco) (3.9) 


where symbol —o denotes a reverse orientation of a spin with respect to o and 
n, is the occupation number of an electron in a quantum state i. 


The occupation numbers are defined now as 


n, = Pole, ayy (E,, EME (3.10) 


with 


1E,.£)=| ef a J- (3.11) 


B 
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being the Fermi Distribution Function (FDF) where E, stands for the Fermi 
level of a crystal; k, is the Boltzmann’s constant and T represents the tem- 
perature in K. 

To find the Fermi energy E, a mean density of states is subsequently de- 
fined 


its 1 
P(E A)= yy PLEA) (3.12) 


with N denoting the number of points of the planar Brillouin zone; y standing 
for the number of electronic bands used in the calculations and also the mean 
number of electrons per atom is given by 


ae 1 
i= 3.13 

“ 2nNy >, na aa 
Now the Fermi level can be calculated as the upper limit of integration in the 
equation 


E; 
7, = | p(E)dE for T>0 (3.14) 


In practice lower and upper limits of integration in (3.10) and (3.14) are both 
finite, chosen in such a way that, outside of this energetic limits, p(E,A)=0 
everywhere. 

The set of occupation numbers n, can be calculated now by means of a self- 
consistent procedure using eqs. (3.6) and (3.10) — (3.14) until the moment when 
a full convergence with assumed accuracy will be achieved in two consecutive 
iterations. 

The quantities n, and E, are fundamental ones for the definition and fur- 
ther calculations of the values reading to the solution of the paradoxes listed 
above. 


4. AN EXPLANATION OF PARADOXES ONTO METALLIC SURFACES 


In order to explain the problems listed in Section 2 it seems to be useful to 
consider still a role of boundary conditions and external polarizations in ECPA. 

In the most cases the metallic sample is “embedded” into an external me- 
dium (vacuum, gas, another metal, electrolyte solution, electric or magnetic 
field, etc.) which acts with electrode surface in self-consistent way. It can be 
generally treated as the boundary condition of the problem. 
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This condition alters the surfacial matrix elements K,,,, in secular equation 
for A, =(I,m,o)=1; A, =(n,m,o)=n by means of adding new element of 
interaction [2, 28] 


Ka Kaa + Wi, = Ki, +W, =K,, tW, (4.1) 


Further in the paper the symmetrical boundary conditions will be used every- 
where, i.e., W, =W, =W. 

The second factor influencing significantly on the LDOS is the polarization 
parameter of a metal caused by and external voltage (denoted as n pı). Polari- 


zation may be defined as the change of a mean value of occupation number 7, 
in electroneutral metal caused by adding to or removing from an atom the elec- 
trons by means of external voltage. It leads, in consequence, to the shift of the 
Fermi level of the crystal. For one-band approximation (only valence-bands 4s, 
5s and 6s for Cu, Ag and Au, respectively) used in this paper, the state 7, =0.5 
stands for an electroneutrality condition (two possible places of electron with 


spins up and down in ns' state), whereas 0<7,,<0.5 and 0.5<n,, <1 denote 


positive and negative charge of the metal, respectively. 

Now we can define in the ECPA terminology two quantities being fundamental 
for the solution of the above mentioned paradoxes of metallic surfaces. These are 
the surface charge of the metal Q, (mentioned earlier in comment to formula (2.1)) 


and corrected work function ¢*, [34]. 


cor 


The surface charge is defined as [1, 2, 29, 34] 


Q = -us| F nn Y-A, uc “cm” (4.2) 


where e stands for the electron charge (without sign); S; represent the number 
of atoms per square centimeter of a surface of the crystallographic face denoted 
as i. Of course, for only valence-band approximation y =1 in eq. (4.2). 

The corrected work function has been derived [1, 34] by means of the ex- 
pansion of the occupation number n, into Taylor’s series. Taking also the 
boundary conditions into account, the final formula for the work function is the 
following: 


03, =-(E, + AE, + E})=—E, - AE, —1,(%, =m) (4.3) 
where the second term of the right-hand side of eq. (4.3) is connected with the 
boundary condition and is always equal to zero for vacuum. The third term in 


(4.3), ie., E’ =1,(7,—n,), describes the surface potential barrier and is di- 
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rectly related to the so-called local work function which can be directly meas- 
ured by means of PAX method [9]. 

Let us now to come back to the list of problems reported in Section 2. 

All the paradoxes can be solved in the framework of GFM and ECPA in the 
following way: 

Ad 1. 

In Table 2 the calculated values of the work function ¢ have been shown 
for Cu, Ag and Au single crystals and compared with the same faces measured 
in experiment, i.e. for (100), (110) and (111) Miller’s faces. The calculations 
for Cu have been also compared with other theoretical methods [11, 33]. It fol- 
lows from a comparison between calculations and averaged values of the meas- 
ured work function from Table 1 that the use of the ECPA gives better agree- 
ment with experiment than EHT (Extended Hiickel Theory) [33] applied by 
Dempsey and Kleinman [35] and much better than Lang and Kohn’s (LK) 
method [11]. The agreement is the best if an average experimental value of the 
work function was obtained from reasonable numbers of different measure- 
ments (five and six cases for Cu(100) and Ag(100) in Table 1, respectively). In 
other cases, particularly when only two different experimental results are avail- 
able such a comparison with experiment is rather poor. 

So generally one can state that GFM and ECPA reproduce quite reasonable 
the corrected work function ¢*, for noble metal monocrystals and one can pre- 
sumably expect that these calculated values are very close to the average ex- 
perimental data for other metals provided that a large enough number of meas- 
urements was carried out. 


Table 2. Corrected work function ¢* calculated for Cu, Ag and Au single crystals of 


(100), (110), (111) faces. In the case of Cu a comparison with other calculations is also 
shown 


¢*,(Cu) [eV] 

| os so tye 100):r ight ood). ir WAEREA Re tee 

ġ2, (Ag) [eV] 
¢2, (Au) [eV] 


a. Ref. [1], b. Ref. [33], i. Ref. [11] 
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Moreover, in all the calculations analysed above, the work function changes 
in the direction f” >! > o (with Me = Cu, Ag, Au) what remains in 
quite good correlation with the most dense surface of (111) face when the high- 
est surface potential barrier occurs in eq. (4.3) and the most open in the case of 
(110) (the smallest barrier). The same tendency appears in Lang-Kohn’s results 
[11] in Table 2 although they remain very far from experimental data. 

In fact, this tendency is only confirmed experimentally for Cu single crys- 
tals. Unfortunately, for Ag and, more particularly, for Au there are big discrep- 
ancies in the experimental data which do not allow us to analyse it. 

Therefore we have shown that the problem of large deviations in the ex- 
perimental data of the work function can be at least partially avoided by the use 
of GFM calculations. But, of course, this conclusion still demands an experi- 
mental confirmation for other metals. 


Ad 2. 

The paradox disappears immediately if we take the corrected work function 
eq. (4.3) instead of the classical one into account. We must put AE, =0 in 
(4.3) because there are no boundary conditions influencing a crystal in a vac- 
uum. Then it is evident that the surface barrier of the potential, i.e. 
AE? = 1, (7, —n,) for a fixed state of a bulk charge (7, =0.5 assuming that the 
measurements are carried out for electroneutrality conditions) depends only on 
n, distribution in the crystal. The value of the Coulomb integral 7, is as- 
sumed to be constant for the nearest neighbours interactions in ECPA [1, 2]. 
For instance, I$" = 1“ =0.4Ry [1, 28]. 

In paper [28] there was calculated a distribution of an occupation number 
along a sample thickness. For Cu(100) and Cu(110) the shape of n, distribution 
is completely different. It means that strong anisotropy of the n, distribution 
takes place also in one sample of a single crystal with different faces. So the 
barrier AE? must be different depending on crystallographic axis of the mono- 
crystals although the Fermi level is only one. 

Ad 3. 

In the light of eq. (4.3) this paradox consists in misunderstanding about an 
assumption p% =0. 

In classical electrochemistry the pzc is understood as a state of compensa- 
tion of a surface charge of an electrode to the zero’s value Q, = 0 (cf. comment 
follows eq. (2.1)). However, it does not mean that a metallic phase lost its elec- 
tric charge in this case. 
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The state with Q, =0 must be treated only as a disappearance of the surface 
potential barrier, so then 7, =n, =n, in eq. (4.3). But such a compensation of 
surface barrier can be only possible by means of an external polarization. 

In the electroneutral state of a metal the surface barrier is equal to 
An=1(0.5—n,) and then the polarization n,, =—An is needed to compensate 


such a barrier. However, it leads to the bulk charge different from electroneu- 
trality, so the metallic phase has to be charged and, in consequence, Ph +0. 


Therefore, eq. (2.1) is correct and the paradox connected with a contact po- 
tential also disappears. 


Ad 4. 

There are only few attempts in electrochemical literature of the last twenty 
years to explain the unusual phenomenon of the ionic specific adsorption. The 
most important theoretical papers [35-39] try to explain this process as a kind 
of chemisorption with covalent bonds between metal and adsorbed ions. These 
bonds seem to be created by means of a charge transfer from a valence band of 
the metal to the partially empty d-bands of halides Cl, Br, T. Needless to say 
that such an explanation takes only a final step of the process into account 
when desolvated ions are already chemisorbed on the metal surface. 

None of the theories gives an answer for the fundamental question of what is 
the reason that two species (electrode surface and ions) with the same positive 
or negative charges attract themselves seemingly against the laws of electro- 
statics. 

An explanation of this question has been provided by means of surface 
charge calculations (eq. (4.2)) as a function of boundary condition W and exter- 
nal polarization 7. 


Table 3, taken from [2], reports the function Q, = f(W,n,,) for Cu, Ag and 
Au single crystals of the face (100). For other cystallographic faces the function 
has similar character [2]. 

There are two regions in Table 3, near the pzc line (from the upper left cor- 
ner to the lower right one), in which the surface and bulk charges of the metal 
reveal the reverse signs. We can see that for n,,=0.3, 0.4 Q,<0 and for 


i, =0.6 0.7 Q, >0. 


From the results analysed above one can conclude that a long-range Cou- 
lombic attraction between the positive bulk of the metal and an anion which is 
coming towards the negative crystal surface determines mainly the nature of the 
specific adsorption forces. Charge transfer is the final step of the process. 
Similar mechanism takes place in the case of the cation adsorption on positively 
charged electrode surface. 
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Table 3. Surface charge 9 [uC/cm?] as a function of w, (in Ry) and 7 E calculated by 
means of ECPA for Ag(100) and Au(100) single crystals. The values of Q, for Cu are 
taken from [1] 


Wi) |a =03 [m,=04 [a =05 [m,=06 [a =07 | 


So an explanation of this very strange phenomenon turned out to be very 
trivial. The superposition of the opposite charges on the surface and in the bulk 
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of the metal produces an electric field attracting ions towards the electrode. 
There is no contradictions in such a behaviour with respect to the classical laws 
of electrostatics. 


In Figure 1, taken from [40], the mechanism of specific adsorption consist- 
ing of three steps is presented . 


= 1 + 
eran 5 te 
1 + j- lns 
by jut - EG — oe 
jae Co i 1 + 
1 + a = ta 
\ { "+ 
q,<0 
Polarizatien>® a 


Te tt StH 


Figure 1. The three-steps mechanism of the specific absorption of anions following from 
the ECPA calculations 


a. Generation of the reverse signed charges in the bulk of metal and on the surface by 
means of an external polarization 


b. Electrostatical attraction of anions by positively charged bulk - adsorption on the 
surface 


c. Electron transfer of valence electrons from the electrode to the empty d orbitals of 
adsorbed anions 


5. CONCLUSIONS 


In this paper the possibilities of the Green Function Method and Extended 
Coherent Potential Approximation of the explanation of some unusual phenom- 
ena in some experiments on metallic surfaces have been reported. 

In practice, this is mainly connected with the work function measurements 
of noble metal single crystals. In this case, when good theoretical results are 
achieved, it is possible to take advantage for reproducing of reasonable tenden- 
cies in the behaviour of the work function as the “vicarious experiment”. 

On the other hand the power of this theoretical method consists in an expla- 
nation of the effects which presumably could be never fully understood only on 
the basis of the experiment regardless of the technical level of the measure- 
ments. It mainly concerns the ionic specific adsorption where an experimental 
separation of the electrode charge between bulk and surface charge of metal is 
practically impossible. The only suitable experimental values are the potential 
of an electrode (measured with respect to certain reference electrode) and its 
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surface charge. As it was shown the GFM and ECPA are able to make such a 
separation in quite a natural way. 

There are also other kinds of controversial effects which may be qualified, 
in fact, as artificial difficulties. For instance, the dependence of the work func- 
tion of a single crystal on the Miller’s indices of its face as well as the paradox 
of contact potential properties. In this case only the good theoretical definition 
of the quantities which are measured (work function) is quite enough to arise an 
agreement between theory and experiment. 

Finally it is worth to emphasise that the GFM’s possibilities are based on the 
level of precision of the Hamiltonian used in the equation of motion (3.3) as well 
as its associated commutation rules of the operators defined in the appropriate 
Fock space. In the bulk the particles or quasi-particles follow only two kinds of 
statistics (bosons or fermions) while on the surface might exist other possibilities 
such as the quantum dots and quantum Hall effect which can lead to the new su- 
perconductors and the giant or even colossal magnetoresistance [41]. 

Very recently, a theoretical paper on the work function calculation by means 
of so-called Brodie’s plasma model has been published [42]. Unfortunately, al- 
though a good agreement of the theoretical results with experimental data (in 
vacuum) for many polycrystallic metals was achieved in this paper there are no 
boundary conditions taken into account. 
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